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Abstract. A detailed radiative transfer analysis of the observed continuum and molecular line emission toward 
the deeply embedded young stellar object IRAS 16293-2422 is performed. Accurate molecular abundances and 
abundance changes with radius are presented. The continuum modelling is used to constrain the temperature 
and density distributions in the envelope, enabling quantitative estimates of various molecular abundances. The 
density structure is well described by a single power-law falling off as r -1,7 , i.e., in the range of values predicted by 
infall models. A detailed analysis of the molecular line emission strengthens the adopted physical model and lends 
further support that parts of the circumstellar surroundings of IRAS 16293-2422 are in a state of collapse. The 
molecular excitation analysis reveals that the emission from some molecular species is well reproduced assuming a 
constant fractional abundance throughout the envelope. The abundances and isotope ratios are generally close to 
typical values found in cold molecular clouds in these cases, and there is a high degree of deuterium fractionation. 
There are, however, a number of notable exceptions. Lines covering a wide range of excitation conditions indicate 
for some molecules, e.g., H2CO, CH3OH, SO, SO2 and OCS, a drastic increase in their abundances in the warm 
and dense inner region of the circumstellar envelope. The location at which this increase occurs is consistent with 
the radius at which ices are expected to thermally evaporate off the grains. In all, there is strong evidence for the 
presence of a 'hot core' close to the protostar, whose physical properties are similar to those detected towards most 
high mass protostars except for a scaling factor. However, the small scale of the hot gas and the infalling nature of 
the envelope lead to very different chemical time scales between low mass and high mass hot cores, such that only 
very rapidly produced second-generation complex molecules can be formed in IRAS 16293-2422. Alternatively, 
the ices may be liberated due to grain-grain collisions in turbulent shear zones where the outflow interacts with the 
envelope. Higher angular resolution observations are needed to pinpoint the origin of the abundance enhancements 
and distinguish these two scenarios. The accurate molecular abundances presented for this low-mass protostar 
serve as a reference for comparison with other objects, in particular circumstellar disks and comets. 

Key words. Stars: formation - Stars: individual: IRAS 16293-2422 - ISM: abundances - circumstellar matter - 
Radiative transfer - Astrochemistry 



1. Introduction 

Low mass protostars in their earliest stages of evolution 
are deeply embedded in large amounts of dust and gas. 
The nature of the emission from such star- forming regions 
makes them ideal to study in the infrared to radio wave- 
length regime. In the last decade, the sensitivity of re- 
ceivers operating at these wavelengths has increased dra- 
matically and the resulting high quality observations, sup- 
plemented by careful modelling, have provided most of the 
current day knowledge about the chemistry and physics, 
and their rapidly changing properties, in early stellar evo- 



lution. While the general scenario of low mass stellar evo- 
lution is reasonably well understood (e.g., Shu et al. 1993; 
Evans 1999; Andre et al. 2000) many uncertainties remain. 
For example, the nature of the hot (T>90 K) and dense 
( n H 2 J>10 7 cm -3 ) regions of gas observed towards some low 
mass protostars is not yet fully established and the possi- 
ble link with the so-called hot cores observed towards most 
high mass protostars needs to be investigated further. 

In the case of high mass star formation, it has be- 
come clear that hot cores represent one of the earli- 
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est phases (Walmsley 1992). Chemically, hot cores are 
characterized by high abundances of fully hydrogenated 
molecules such as water (H2O), ammonia (NH3) and hy- 
drogen sulfide (H 2 S), along with a rich variety of com- 
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plex organic molecules ranging from methanol (CH3OH) vey presented in Blake ct al. (1994) and van Dishocck 



and et jianol (CH 3 CH 2 OH) to methyl cyanide (CH 3 CN), 
dimethyl ether (CH3OCH3), m ethyl formate (HCOOCH 3 ) 
and ethyl cyanide (C 2 H 5 CN) (|Wahnsley fc Schilke 1993| ; 
Kuan fc Snyder 1996| ; |Hatchell et al. 1998] ; |Schilke 2000| ). 
The chemical richness is explained by evaporation of the 
ice mantles above ~90 K, followed by rapid gas-phase 
ion-molecule reactions leading to more complex species 
for a period of >10 4 yrs (Charnley et al. 1995; Millar et 
al. 1997; Rodgers & Charnley 2001, see van Dishoeck & 
Blake 1998 and Langer et al. 2000 for overviews). Low 
mass protostars are less luminous and less massive, but 
a simil ar physical structure is expected except for a scale 
factor flCeccarelli et al. 1996| ; [ivezic fc Elitzur 1997| ). On 
the other hand, shocks due to the interaction of the out- 
flows with the envelope can also liberate ice mantles and 
drive a high-temperature chemistry; such shocks may be 
relatively more important for low-luminosity objects than 
for high- mass protostars (e.g., van Dishoeck et al. 1995). 
It is of considerable interest to establish if low mass proto- 
stars also have hot and dense regions, and if so, whether a 
similarly complex organic chemistry to that found in the 
case of high-mass protostars has ensued and whether pas- 
sive heating by the accretion luminosity or active shocks 
dominate the liberation of grain mantles. Since it is the 
material in the warm inner envelope that will be incor- 
porated into circumstellar disks, it is important to know 
the level of chemical complexity as it relates to forming 
planetary systems. 

IRAS 16293-2422 is by far the best candidate for in- 
vestigating a low mass hot core (e.g., Blake et al. 1994; 
van Dishoeck et al. 1995; Ceccarelli et al. 2000a,b). IRAS 
16293-2422 is a deeply embedded low mass protostellar 
object located within the p Ophiuchus molecular cloud 
complex. Due to the relative proximity of this source 
(160 pc; Whittet 1974 ) a wealth of molecular lines has 
been de tected, in spite of i ts relatively low luminosity 
(27 L ; |Mundy et al. 1986j ), and this has made IRAS 



et al. (1995). It was found that molecular line emission is 



16293-2422 one of the best studied young stellar objects 
Interferometer observations of radio and millimetre con- 
tinuum emission reveal two c ompact source s in the center 
of it s circumstellar envelope ( Wootten 1989 ; Mundy et al 
1990, 1 992| ; Looncy ct al. 200C ), likely to be accretion disks 
through which matter is fed onto the central stars. The 
s eparation of the tw o protostars is approximately 800 AU 
([Looncy et al. 2000|). IRAS 16293-2422 is thought to be in 



potentially a very powerful tool to probe both the physics 
and chemistry of the circumstellar environment; however, 
a full radiative transfer analysis was not performed and the 
derived abundances have significant uncertainties. Due to 
the complexity of molecular excitation and its sensitivity 
to the environment, various species — and even different 
lines of the same species — probe different parts of the cir- 
cumstellar material. At least three physically and chem- 
ically distinct parts were identified including a circumbi- 
nary envelope, circumstellar disk(s), and outflow compo- 
nents. The latter component was thought to be a small 
and warm region of a few arcsec in size where the bipolar 
outflow(s) interact with the inner part of the circumbinary 
envelope. Recently, Ceccarelli et al. (2000a, b) used deep 
JCMT observations of H2CO combined with a physical- 
chemical collapse model to argue that IRAS 16293-2422 
does in fact have a hot-core-like region in which the lib- 
eration of ices is consistent with heating by the accretion 
luminosity. 

We present here spherically-symmetric radiative trans- 
fer modelling of the dust and gas components constitut- 
ing the material in the circumstellar envelope of IRAS 
16293-2422. The dust parameters are constrained by the 
observed continuum emission in the form of submillime- 
tre brightness maps and the SED over a large wavelength 
region. The resulting temperature and density structures 
are a prerequisite to chemical studies of the molecular gas 
present in the envelope. The approach taken here is dif- 
ferent from that adopted by Ceccarelli et al. (2000a,b) in 
that the physical parameters of the envelope are derived 
empirically from the analysis of the dust emission. Once 
the physical structure of the envelope is known, a detailed 
excitation analysis of molecular millimetre line emission is 
performed aimed at obtaining accurate abundances. This 
provides valuable insight into the complex chemistry oc- 
curring in this proto-stellar envelope. In particular, the de- 
rived abundances allow for direct comparison with other 
sources and comets (e.g., Bockelee-Morvan et al. 2000). 
Moreover, searches for evidence of abundance changes, 
e.g., due to evaporation of ices ('jump models') is of con- 
siderable interest. In addition to providing constraints on 
the chemistry, the molecular line observations give further 
information on the physical structure, e.g., kinematic in- 



formation. Similar strategies have been adopted by van 



one of :he earliest stages of formation; the observed spec- 



tral en ;rgy distribution (SED) can be fitted by a modified 



blackbody of ~40K (e.g., |Walker et al. 1986| ; [Andre et al 
2000|)~i|gid has a high ratio of submillimetre to bolomctric 
luminosity, suggesting a large amount of envelope mass. 
This places IRAS 16293-2422 among a family of deeply 
embedded and recently formed hydrostatic stellar objects 
known as, in the traditional evolutionary sequence of low 
mass protostellar objects (e.g., Andre et al. 1993), 'class 
0' protostars. 

The circumstellar surroundings of this protobinary 
star were extensively studied in a large molecular line sur- 



dcr Tak et al. (1999| , |2000b| ) and [Hogcrheijde fc Sandcf 
(20001 ) and have proven to be powerful tools when deter 



mining the physics and chemistry of star-forming regions. 

2. Observations and data reduction 

In this Section the observational constraints used in the 
radiative transfer modelling of the circumstellar envelope 
around IRAS 16293-2422 are presented. The SED and 
submillimetre continuum brightness maps constrain the 
physical structure of the envelope. Millimetre molecular 
line observations provide further information on the phys- 
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Table 1. The spectral energy distribution of IRAS 16293- 
2422. 



A -Fobs' 

[/im] [Jy] 



AF ol 



-frnod 



[Jy] ["] Ref. [Jy] 



2900 


0.60 


0.13 


60.0 


3 


0.58 


1300 


6.97 


2.24 


30.0 


3 


7.5 


850 


30.8 


6.2 


15.0 


1 


31.3 


450 


220.4 


44.5 


8.8 


1 


220.3 


100 


1032.0 


226.0 


237.0 


2 


852.4 


60 


254.9 


59.5 


160.0 


2 


162.4 



a Observed flux integrated over emitting region. 
b A calibration uncertainty of 20% has been added 
c Size of the main beam. 

d Flux predicted by best fit model using a density structure 
described by a single power-law and OH5 dust opacities. Note 
that the 60 /im and 2.9 mm fluxes are not included in the mod- 
elling (see text for details). 

Refs. - (1) This pa per; (2) IRAS Point Source Catalogue; (3) 
Walker et al. (1990h. 
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Fig. 1. SCUBA images at 450 and 850 //m of IRAS 16293- 
2422. The contours start at the 3cr level (0.9 and 0.24 
Jybeam -1 for 450 and 850 /im, respectively) and increase 
by multiples of 2. 



ical structure, in particular the large scale velocity field, 
and allow for studies of the chemistry in the envelope. 

2.1. Spectral energy distribution 

The SED of IRAS 16293-2422, as presented in Table [l] is 
used to constrain the total amount of material present in 
the circumstellar envelope. Care has been taken that only 
total fluxes integrated over the whole emitting region for 
wavelengths <1.3 millimetre (mm) are used in the SED 
modelling to ensure that emission from the envelope itself 
is dominating the observed flux. The contribution from 
the disk(s) to the observed emission starts to become im- 
portant at wavelengths longer than ~lmm. However, any 
disk emission would be significantly diluted in the rela- 
tively large beams used here (e.g., Motte & Andre 2001). 
That thermal emission from an extended dusty envelope is 



responsible for the majority of the observed millimetre flux 
is confirmed by the radiative transfer analysis performed 
in Sect. ||. The model fluxes predicted by the best fit en- 
velope model are compared with observations in Tabic [l| 
Additionally, the observed IRAS flux at 60 /im is not used 
in the analysis since the emission at that wavelength em- 
anates from the inner hot parts of the envelopes where 
the dust grain properties are probably significantly differ- 
ent from those in the cooler outer parts. Dust grains in 
the outer parts will be coated by a layer of ice which, as 
the temperature increases towards the star, starts to evap- 
orate thereby changing the optical properties of the dust 
grains. In Sect. |I| the effect of varying the dust opacities 
will be investigated. A relative calibration uncertainty of 
20% is added to all flux measurements, in addition to any 
statistical errors, which dominates the error budget and 
gives all points on the SED more or less equal weights. 

2.2. Submillimetre continuum observations 

Submillimetre continuum observations were retrieved 
from the James Clerk Maxwell Telescope (JCMT) public 
archive^. The data were obtained during an observational 
run in April 1998, using the Submillimetre Common-User 
Bolometer Array (SCUBA), and consist of pairs of images 
at 450 /xm and 850 /im. The dual SCUBA array contains 
91 pixels in the short-wavelength array and 37 pixels in 
the long-wavelength array, each covering a hexagonal 2.'3 
field. The SC UBA bolometer array camera is described in 
some detail in Holland et al. (1999 ). 

The imaging was made using the jiggle-mapping mode 
to produce fully sampled maps. In this mode the SCUBA 
bolometers instantaneously under-sample the sky and a 
64 point jiggle pattern is carried out by the telescope to 
fully sample both the long and short wavelength arrays. 
In practice, the 64 point pattern is broken down into four 
16-point sub-patterns that spend Is integrating on the 
source at each gridpoint. After a sub-pattern has been 
completed the telescope is nodded and the pattern is re- 
peated again so that sky subtraction can be made. In all, it 
takes 128 s to complete one full on/off source jiggle map. 
The jiggle mapping mode is the preferred observational 
mode for SCUBA when imaging sources smaller than the 
chop throw. Usually, a chop throw of 2' is used to ensure 
chopping off the arrays. A larger chop throw would result 
in poor sky subtraction and loss in image quality. 

The data were reduced in a standard way, as de- 
scribed in Sandell (1997), using the SCU BA reduction 



package SURF ( Jenness fc Lightfoot 1997 ). The images 
were calibrated using simultaneous observations of Uranus 
retrieved from the JCMT archive. The sky opacities at 



tittp : //www. jach.hawaii . edu/ JACpublic/ JCMT/ 



The JCMT is operated by the Joint Astronomy Centre in Hilo, 
Hawaii on behalf of the present organizations: the Particle 
Physics and Astronomy Research Council in the United 
Kingdom, the National Research Council of Canada and the 
Netherlands Organization for Scientific Research. 
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Table 3. Additional molecular line observations of IRAS 16293-2422 using the JCMT. 



Frequency / T mb dV a T mb b AV b 



Molecule 




Line 


[MHz] 


[K kms" 1 ] 


[K] 


[kms 


- 1 ] Ref. 


13 CO 


J 


= 2-1 


220401.7 


50.0 


9.7 


3.5 


1 




./ 


= 3-2 


330588.1 


67.2 


8.5 


4.7 


1 


c 18 o 


J 


= 2-1 


219560.4 


21.6 


7.8 


2.6 


1 




J 


= 3-2 


329335.0 


33.6 


10.1 


3.2 


1 


c 17 o 


./ 


= 2-1 


224714.4 


6.6 


2.5 


2.7 


1 




J 


= 3-2 


337061.1 


10.9 


3.0 


3.4 


1 


cs 


./ 


= 5-4 


244935.7 


41.1 


11.1 


3.3 


1 




./ 


= 7-6 


342887.8 


51.4 


13.6 


3.3 


1 


c 34 s 


./ 


= 5-4 


241016.2 


4.8 


1.3 


3.4 


1 




./ 


= 7-6 


342883.0 


5.7 


1.6 


3.2 


1 


HCN 


J 


= 3-2 


265886.4 


51.3 


4.0 


7.5 


1 




J 


= 4-3 


354505.5 


63.4 


5.0 


7.3 


1 


HNC 


J 


= 3-2 


271981.1 


14.2 


4.7 


2.9 


1 




J 


= 4-3 


362630.1 


11.9 


4.1 


3.0 


1 


HCO+ 


J 


= 3-2 


267557.6 


70.3 


18.0 


3.5 


1 




J 


= 4-3 


356734.0 


95.9 


21.5 


4.0 


1 


H 13 CO+ 


J 


= 3-2 


260255.5 


10.3 


3.8 


2.6 


2 




J 


= 4-3 


346998.5 


8.5 


3.1 


2.7 


1 


29 SiO 


J 


= 8-7 


342979.1 


1.8 


0.3 


5.4 


1 


HC 3 N 


J 


= 24-23 


218324.8 


0.64 


0.15 3.9 


1 



a Total integrated intensity calculated over full extent of line. The calibration uncertainty in the intensity scale is estimated to 
be ~l5-20%. 

b Estimated from a gaussian fit to the observed spectrum. A colon (:) indicates an uncertain value due to low signal-to-noise 
or, in the majority of cases, a significant departure from a gaussian line profile. 
Refs. - (1) JCMT public archive; (2) This paper. 



Table 2. Observational results from SCUBA images. 



A 


r a F tot b 


-fpeak 


rms 


#mb d 


H 


[Jy] 


[Jybeam - ] 


[Jybeam - ] 


["] 


450 


0.5 - 0.7 220.4 


72.0 


0.3 


8.8 


850 


0.12 - 0.15 30.8 


16.8 


0.08 


15.0 



a Zenith opacity of the atmosphere. 

b Total flux integrated over emitting region. The calibration 
uncertainty is estimated to be ~20%. 

c Flux at stellar position. The calibration uncertainty is 
estimated to be ~20%. 



Geometrical mean of the beam size. 



450 /mr and 850 /im were estimated using the 1.3 mm opac- 
ity, monitored by the Caltech Submillimcter Observatory 
and listed for each individual observation, using the re- 
lations in Archibald ct al. (200C). The validity of these 



relations have been checked and confirmed using SCUBA 
sky dips. The total calibration uncertainty is estimated to 
be approximately ±20% at 450 /im and about ±10% at 
850 /im. Care was taken to select data taken during good 



mined from the Uranus observations and is dominated by 
an approximately gaussian main beam with a deconvolved 
FWHM of 14'.'2 x 16"0 and 8'.'5 x 9"1 at 850 /im and 450 /mr, 
respectively. A substantial error beam is, however, present 
at these wavelengths (see below) picking up significant 
amounts of flux. The error lobe pick up is estimated to 
be approximately 15% and 45% at 850 /im and 450 /im, 
respectively, and is taken explicitly into account in the 
analysis. 

The final 450 /im and 850 /im images are presented 
in Fig. [l]. All offsets reported are relative to the 
adopted position of the protobinary star IRAS 16293- 
2422 (cv 20 oo = 16 h 32 m 22!91, <5 2 ooo=-24°28'35"6). The lo- 
cations of the two protostars IRAS 16293A (MM1) and 
IRAS 16293B (MM2) relative to this position are (-3", 
— 1") and (—5", ±3"), respectively. The emission appears 
to have an overall spherical symmetry and is centered on 
the adopted central position within the pointing accuracy 
of the telescope. For comparison, the pointing accuracy of 
the JCMT is estimated to be about ±1.5" in both eleva- 
tion and azimuth. Also visible is a second, weak, compo- 
nent apparent near the eastern edge of the SCUBA maps. 
This second component, or IRAS 16293E (+77", -22"), 
which was first identified by its strong ammonia emission, 



is most probably also a class protostar (Mizuno et al 



to excellent submillimctrc conditions. The beam is deter- 1990; Castets et al. 2001) 
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The FWHM of the emission centered on the stel- 
lar position is 20'.'8x 19'.' 4 and 21'.'9x 19'.'3 at 850 ^iia and 
450 /im, respectively. The deconvolved envelope sizes as- 
suming both the beam and brightness distribution to be 
described by gaussian functions are ^14" at 850 ^m and 
~19" at 450 /im. Thus, only the 450 /jm emission appears 
to be resolved. The observational results are summarized 
m Table |. To compare the observed brightness distri- 
butions with the predictions from a spherically symmet- 
ric model, the SCUBA maps were azimuthally averaged 
in bins with half the corresponding beam size in width. 
Moreover, care was taken to block out any contribution 
from IRAS 16293E. The resulting radial brightness distri- 
butions are shown in Fig. ||. 

2.3. Millimetre molecular line observations 

A survey of the millimetre molecular line emission to- 



wards IRAS 16293-2422 was presented in Blake et al 
|(1994p and |van Dishoeck et al. (1995| ). This large data 



set forms the base for the molecular excitation analysis 
performed in this paper. The absolute calibration uncer- 
tainty of the intensities is estimated to be ~30%. In addi- 
tion, we have searched the JCMT public archive for com- 
plementary millimetre line observations. This additional 
set of data, taken at face value, is presented in Table [| 
Lines for which multi-epoch observations are available in 
the JCMT archive typically display intensities that are 
co nsistent to ^20%. This was also the conclusion reached 
by |Schoier fc Olofsson (200l] ) for a large survey of car- 
bon stars. When newer data were available for a particu- 
lar transition, they were usually adopted, due to a higher 
signal-to-noise and/or greater spectral resolution. In one 
case [H 13 CO + ( J=3^2)] the old data set was found to 
have a significantly lower line intensity, possibly due to 
pointing problems. Additional H2CO data published re- 
cently by Loinard et al. (20"(X| ) were further used in order 
to increase the number of observed transitions for this 
molecule. 

The detected molecular line emission probes the full ra- 
dial range of the envelope, providing additional constraints 
on the physical structure of IRAS 16293-2422. Only infor- 
mation on the lowest transitions of the molecules, which 
occur at millimetre wavelengths and probe the very cold- 
est outer parts, is lacking. The observed line shapes pro- 
vide valuable information on the velocity structure in 
the envelope. For example, the single-dish observations of 
abundant molecules like CO and CS typically show lines 
with strong self-absorption and some degree of asymme- 
try. The variation of the line profiles among the CO and 
CS isotopomers is potentially a sensitive prob e of infall 
models and will be further investigated in Sect. 4.3. 



3. Radiative transfer 

In this Section the continuum and molecular line radiative 
transfer codes used to constrain the physical and chemi- 
cal structure of the envelope around IRAS 16293-2422 



are presented. The envelope is assumed to be spherically 
symmetric and the approach adopted here is to first de- 
termine the density and temperature structures from the 
dust modelling. This, in turn, allows for abundances of 
various molecules present in the circumstellar envelope to 
be determined. 

3.1. Dust radiative transfer model 



In order to model the observed continuum emission and to 
be able to extract some basic parameters of the dusty en- 
velope around IRAS 16293-2422 the p ublically available 
dust radiative transfer code DUSTY^ fllvczic et al. 1999| ) 
has been adopted. DUSTY makes use of the fact that 
in some very general circumstances the radiative transfer 
problem of th e dust possesses scaling properties ( tvezic fe 
Elitzur 1997). The solution is presented in terms of the 
distance, r/n, scaled with respect to the inner boundary 
n. In addition to the properties of the dust and the rel- 
ative size of the envelope, the only parameter needed for 
a full description of the problem is the spectral shape of 
the radiation emitted by the central source. This means 
that the luminosity is totally decoupled from the radiative 
transfer problem and it is only used to scale the solution in 
order to obtain the absolute distance scale. This scaling 
property of the dust radiative transfer is very practical, 
in particular when modelling a large number of sim ilar 
sources, and is put to use in forgensen et al. (2002 ) for 
a survey of low mass protostars at various evolutionary 
stages and in Hatchcll et al. (2000 ) to model a selection 
of high mass protostars. 

The most important parameter controlling the output 
is the dust optical depth 



t\ = K\ / pd{r) dr, 



(1) 



where k\ is the dust opacity, pa the density distribution of 
the dust in the envelope, and n and r e the inner and outer 
radius of the dust envelope, respectively. Introducing the 
dust-to-gas mass ratio, S, allows Eq. [I] to be written 



n H2 (r) dr = K A £>m H2 A ? H 2 



(2) 



where nn 2 is the number density distribution of molecular 
hydrogen, mu 2 the mass of an hydrogen molecule, and A^h 2 
the column density of molecular hydrogen. In the deriva- 
tion of Eq. I any possible drift velocity between the dust 
and the gas has been neglected. In what follows 5=0.01 is 
assumed. The dust opacities from Osscnkopf fc Hcmiing] 



(1994) were used, corresponding to coagulated dust grains 
with thin ice mantles, at a density of riH 2 ~10 6 cm -3 (col- 
umn 5 in their Table 1, hereafter OH5). Van der Tak et al. 
(1999) considered various sets of dust optical properties 
when modelling the high mass young stellar object GL 
2591, and found that models using the OH5 opacities were 
the only ones that gave envelope masses consistent with 

2 http : / / www . pa . uky . edu/ "moshe/dusty/ 
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those derived from the modelling of the molecular line 
emission. In Sect. 4.2 another set of opacities will also be 
considered, more appropriate for regions where the ices 
have evaporated from the dust grains. 

The dust temperature at the inner radius of the en- 
velope is fixed to 300 K and this sets the inner radius. 
The choice of this temperature is motivated by the ob- 
servations of line emission arising from highly excited 



molecules in the envelope. However, the inner regions 



are complex with breakdown of spherical symmetry and 
interactions between disk, envelope and outflow. In the 



3.2. Molecular line radiative transfer model 

In order to derive accurate molecular abundances for the 
wealth of molecular line emission detected toward this 
source the de tailed non-LTE radiative transfer code of 
[Schoier (2000 ), based on the Monte Carlo method, was 
used. The code produces output that is in excellent agree- 
ment with the Monte Carlo code presented by Hogerheijde 
fc van der Tak (2000] ). It has also been tested against other 



molecular line radiative transfer codes, for a number of 



benchmark problems, to a high accuracy (van Zadelhofi 



preset analysis these complications are ignored and for ct al. 2002) 



simplicity a smooth and spherically-symme tric envelope 
is assumed. Recently, Ceccarelli et al. (2000a ) successfully 
modelled molecular line emission in the envelope around 
IRAS 16293-2422 down to -30 AU, assuming spherical 
symmetry. While the se paration of the two protostars is 
approximately 800 AU ( Looney et al. 2000 ), one of the 
protostars IRAS 16293A (MM1) exhibits jet-like centime- 
tre wavelength emission, water maser emission and associ- 
ated millimetre molecular emission, thus appearing to be 
significantly more active (Wootten & Loren 1987; Mundy 
et al. 1992; Schoier et al. 2002a, in prep.), which justifies 
the approach taken here. The central source of radiation 
is assumed to arise from a blackbody at 5000 K. This is an 
oversimplification considering the binary nature of IRAS 
16293-2422 and the uncertainty in the intrinsic SED of a 
protostar. The final model does not depend on the exact 
stellar temperature adopted, however, within a reasonable 
range of values, since the radiation is totally reprocessed 
by the circumstellar dust. The input parameters are sum- 
marized in Table [|. 

The observational constraints, as presented in Sect, g, 
are the SED and radial brightness distributions at 450 /im 
and, to a lesser extent, 850 /im. The ability of the model 
to reproduce the observational constraints are quantified 
using the chi-squared statistic 



N r 

* 2 = E 



(Ft 



mod.i 



-^obs,i) 



(3) 



where f is the flux and Oj the uncertainty in observation i, 
and the summation is done over all N independent obser- 
vations. The radial brightness distributions from DUSTY 
are extended into 2D surface brightness maps and con- 
volved with the beam as determined from planet obser- 
vations. The beam convolved maps are then azimuthally 
averaged in 3" bins. In the x 2 fitting procedure only data 
points separated by one full beam are used since the \ 2 ~ 
analysis, in practice, requires uncorrelated measurements. 
In the analysis only the inner 50" of the brightness dis- 
tributions are considered which are well above the back- 
ground emission and should not be significantly affected 
by the 120" chop throw. Furthermore, the IRAS fluxes 
from the model were convolved with the proper filters be- 
fore the x 2 analysis of the SED was made. The results 
from the dust radiative transfer are presented in Sect. [|. 



Adopting the parameters of the circumstellar enve- 
lope derived from the dust radiative transfer analysis, the 
Monte Carlo code calculates the steady-state level popu- 
lations of the molecule under study, using the statistical 
equilibrium equations. In the Monte Carlo method, infor- 
mation on the radiation field is obtained by simulating 
the line photons using a number of model photons, each 
representing a large number of real photons from all tran- 
sitions considered. These model photons, emitted locally 
in the gas as well as injected from the boundaries of the 
envelope, are followed through the envelope and the num- 
ber of absorptions are calculated and stored. Photons are 
spontaneously emitted in the gas with complete angular 
and frequency redistribution, i.e., the local emission is as- 
sumed to be isotropic and the scatterings are assumed to 
be incoherent. The weight of a model photon is continu- 
ously modified as it travels through the envelope, to take 
the absorptions and stimulated emissions into account. 
When all model photons are absorbed in, or have escaped 
from, the envelope the statistical equilibrium equations 
are solved and the whole process is then repeated until 
some criterion for convergence is fulfilled. Once the molec- 
ular excitation, i.e., the level populations, is obtained the 
radiative transfer equation can be solved exactly. The re- 
sulting brightness distribution is then convolved with the 
appropriate beam to allow a direct comparison with ob- 
servations. In this analysis, the kinetic te mperature of the 



gas is assumed to follow that of the dust (Ceccarelli ct al 



1996| ; [Doty fc Ncufeld 199J ; |Ceccarelli ct al. 2000a| ) 



Typically, energy levels up to ~500K in the ground vi- 
brational state were retained in the analysis. Vibrationally 
excited levels are not included since the radiative excita- 
tion due to the dust is generally inefficient and, for the 
temperature and density ranges present here, collisional 
excitation to these levels is negligible. Line emission from 
rotational transitions within vibrationally excited states 
has, however, been observed for some species [e.g., the 
GS(u=l, J=7^6) ([Blake et al. 1994j )]. In Sect. | the na- 
ture of such emission is discussed further. Collisional rate 
coefficients are taken from the literature and, in the case of 
some linear molecules, extrapolated both in temperature 
and to transitions involving energy levels with higher J 
quantum numbers when needed [see Schoier et al. (2002b, 
in prep.) for details]. For other (non-linear) species for 
which such extrapolations are not obvious, the excitation 
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is assumed to be in LTE for levels for which no collisional 
rate coefficients are available. 

In what follows, all quoted abundances, fx, for a par- 
ticular molecular species X are relative to that of molecu- 
lar hydrogen, i.e., 



Mr) 



nx(r) 



(4) 



and are initially assumed to be constant throughout the 
envelope. Subsequently, it will be shown that in order to 
model the observed line emission for some molecules, e.g, 
H 2 CO and SiO, this latter constraint has to be relaxed and 
a jump in / needs to be introduced (see also Ceccarelli et 
al. 2000a,b). Van der Tak et al. (2000b) used a similar 
approach for the case of massive protostars. 

The beam profile used in the convolution of the mod- 
elled emission is assumed to be gaussian which is appro- 
priate at the frequencies used here. The best fit model 
is estimated from the x 2 -statistic defined in Eq. || using 
the observed integrated intensities and assuming a 30% 
calibration uncertainty. Even though the new data pre- 
sented in Table || appear to be slightly better calibrated, 
the old data set constitutes the vast majority of obser- 
vational constraints so the larger calibration uncertainty 
is applied to the full data set, for simplicity. Only in the 
cases of CO and CS was a lower calibration uncertainty 
of 15% adopted. These molecules are regularly observed 
towards IRAS 16293-2422 and used as standard spectra. 

4. The physical structure of the envelope 

The envelope parameters, such as the density and tem- 
perature structures, are constrained in this Section using 
mainly the observed continuum emission. 



4.1. Simple analysis 

For optically thin dust emission one can derive the fol- 
lowing expression for its intensity, in the Rayleigh- Jeans 
limit, as a function of the impact parameter b (Shirley 
|et al. 2000[ ) 



= — ,7 = « + 0-l, 



^(0) 

assuming that the density 

and temperature 
T d (r) oc r~ & 



(5) 



(6) 



(7) 



follow power-law distributions. For IRAS 16293-2422 
7~1.8 is derived in the range 6~15" — 50" for the 450 /im 
data. Assuming ,5=0.4 (c.f., Doty & Leung 1994; Shirley 
et al. 2000; see also Fig. |^), a^2.4 is obtained, consis- 
tent with the range of values determined by Shirley et 
al. for a sample of protostellar objects. As discussed by 
Shirley et al. there are a number of caveats when using 



this simplistic approach, e.g., the validity of applying the 
Rayleigh- Jeans approximation in the cool outer parts of 
the envelope and deviations of the dust temperature from 
a single power-law (see also Hogerheijde & Sandell 2000; 
Doty & Palotti 2002). Shirley et al. estimate the uncer- 
tainties in the derived a to be of the order ±0.5 using this 
simple approach. 

4.2. Power-law density model 

In order to derive more reliable envelope parameters the 
dust radiative transfer model presented in Sect. 3.1 is used. 
The density of the dust (and gas) is assumed to follow a 
simple power-law 



1000 AU 



(8) 



where uq is the number density of H2 at a distance of 
1000 AU from the star. The models where the density 
structure is given by a single power-law will be referred 
to as static envelopes, since in the molecular excitation 
analysis no large-scale velocity field is included. The total 
amount of dust and its spatial distribution are simultane- 
ously determined, i.e., T100 (the optical depth at 100 /im), 
a, and r e /n (the geometrical envelope thickness) are the 
adjustable parameters in our model. Scaling the model 
output to the source luminosity of 27 Lq and a distance 
of 160 pc fixes the absolute scale and allows various phys- 
ical parameters of the envelope around IRAS 16293-2422 
to be determined. 

To find the best fit model in the 3D parameter space 
the strategy adopted is to first analyze the model grid by 
applying the observed radial brightness distributions. As 
shown in Fig. |^, the 450 urn and 850 /im brightness distri- 
butions are sensitive to the slope of the density distribu- 
tion and, to a lesser extent, the size of the dusty envelope. 
Changing the total amount of dust by changing tioo has 
only a minor effect on the allowed values of a since nor- 
malized radial brightness distributions are used. Values of 
a in the range 1.5—1.9 are found to be acceptable, with a 
preferred value of 1.7. A typical acc uracy of ±0.2 in the de- 
rived value of a was also found by Jgrgensen et al. (2002 ) 
when analyzing SCUBA images at 450 /im and 850 /im for 
a large sample of protostars. In general, the 450 /im data 
should provide a more reliable value of a because of the 
higher resolution, which provides a larger sensitivity to 
changes in the density structure. 

The SED provides a good constraint on r once the 
density profile is known (Fig. ||; see also see also Doty 
& Palotti 2002). For a density slope of 1.7 the optical 
depth at 100 /im is estimated to be approximately 4.5. At 
450 /im and 850 /im the optical depths are ~0.4 and ~0.1 
respectively. For a>1.7 the size of the envelope is not well 
constrained which is not surprising given that only points 
in the brightness distributions out to 50" (r c /rj~250) are 
used. The circumstellar envelope will eventually merge 
with the more extended cloud material in which the object 
is embedded. The maximum outer radius of the envelope 



8 



F. L. Schoier et al.: Does IRAS 16293-2422 have a hot core? 



Table 4. Summary of the dust radiative transfer analy- 
sis of IRAS 16293-2422 using a single power-law density 
distribution (see text for details). 



Fixed input parameters 

Distance, d 160 pc 

Luminosity, L 27 Lq 

Stellar temperature, T* 5000 K 

Dust temperature at n, Td(n) 300 K 

Dust opacity (OH5) at 100 /im, kioo 86.5cm 2 g _1 

Variable input parameters 



Dust optical depth at 100 /J,m, tioo 
Density power law index, a 
Envelope thickness, r c /ri 



3.0-6.0 
1.4 - 2.0 
50 - 300 



Best fit parameters 



Dust optical depth at 100 /jm, tioo 
Density power law index, a 
Envelope thickness, r c /n 



4.5 
1.7 

250 



Derived parameters 



Inner envelope radius, n 


4.8 x 10 14 cm 


Outer envelope radius, r c 


1.2 x 10 17 cm 


H2 column density, NQI2) 


1.6 x 10 24 cm" 


-2 


H 2 density at 1000 AU, n 


6.7 x 10 6 cm" 


-3 


Envelope mass, M env 


5.4M 




Bolometric flux, fboi 


3.4 x 10' u Wm" 





is fixed at the point where the dust temperature reaches 
10 K. The envelope size r e /n is estimated to be 250 for 
the adopted a of 1.7. 

The quality of the best fit model can be judged from 
the reduced x 2 obtained from 



Xrcd 



y 2 . 
A mi 11 



(9) 



N -p 1 

where p is the number of adjustable parameters. It is found 




850^m 
T, nn =4.5 




Fig. 2. x 2 maps showing the sensitivity of the single 
power-law model to the adjustable parameters T100, ct, and 
Ai?=r c /n using the SED and the SCUBA maps as obser- 
vational constraints. Contours are drawn at x 2 nin +(2.3, 
4.6, 6.2) indicating the 68% ('lcr'), 90%, and 95% ('2<r') 
confidence levels, respectively. 

same radius, in excellent agreement. Thus, IRAS 16293- 
2422 appears to have one of the most massive envelopes of 



the known class protostars (Andre et al. 200C; l0rgensen 



that th |e best fit model presented above has xt^A— 0-1 when 



et al. 2002) 



the SED is used as constraint and Xr O d~0-4 f° r ^ ne 450 M m 
and 850 /im radial brightness distributions, respectively. 
In all, the observations are well reproduced. The best fit 
solution is presented in Table ^. The inner radius, n, of 
the envelope is calculated to be 4.8xl0 14 cm (32 AU) fix- 
ing the outer radiu s to 1.2xl0 17 cm (8000AU or 50"). 
Mundy et al. (199C) estimated the radial size of the enve- 
lope to be ~4400 AU from observations of ammonia emis- 
sion. The best fit model is presented in Fig. [3] overplotted 
on the observational constraints. The total mass of molec- 
ular hydrogen contained within the outer radius of 50" 



Although disk emission is typically only responsible for 
a small fraction of the total flux at sub-millimetre wave- 
lengths it can contribute to the fluxes of the innermost 
points on the brightness profiles leading to a steeper in- 
ferred density profile. Tests in which the flux within a 
radius of one beam was reduced by 50% indicate that the 
best fit value of a is reduced by 0.1-0.2. 

The temperature and density structures obtained from 
the best fit model are presented in Fig. |^. For comparison, 
the predicted temperature structure based upon an opti- 



is estimated to be 5.4 Mq. Blake et al. (1994) derived a 
H2 mass content of ^0.5— 0.75 Mq within a 10" radius 
from an excitation analysis of the observed C 17 molecu- 
lar line emission assuming its abundance to be 3.8xl0~ 8 . 
The present analysis gives a mass of ~0.7Mq within the 



cally thin approximation (Chandler & Richer 200C) and 
scaled to the luminosity of IRAS 16293-2422 is shown for 
two different opacity laws. Both predict the dust temper- 
ature to follow a single power-law. Clearly, the tempera- 
ture structure obtained from the detailed radiative trans- 
fer analysis is not well described by a power-law and has 
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Fig. 3. Best fit model, using a density structure described by a single power-law, compared with observed radial 
brightness distributions and the SED. In the x 2_ana ly s i s of the radial brightness distributions only the data points 
separated by one full beam width were used, shown here with error bars. The error bars represent the rms scatter of 
the observations within each of the bins and are a combination of the noise as well as gradients and departure from 
spherical symmetry in the brightness map. The dotted line is the azimuthally averaged SCUBA beam at the time of 
the observation. In the SED panel the observations, represented by circles with error bars (open circles were not used 
in the x 2 -analysis) , are overlayed with the output from the best fit model. 



a significantly steeper gradient in the inner parts of the 
envelope where optical depth effects are important. The 
interstellar radiation field is potentially important for the 
temperature structure in the outer parts of the envelope. 
However, detailed modelling (S. Doty, priv. comm.) shows, 
for the envelope around IRAS 16293-2422, that this ef- 
fect is small when assuming a typical interstellar radiation 
field. The difference in the dust temperature is ~10-20% 
(a few K) at r>8xl0 16 cm. 

The dust opacities adopted in the previous analysis 
might not be appropriate in the inner hot parts of the 
envelope where the ices start to evaporate off the grains. 
Instead opacities for bare grains, without any ice mantles, 
should preferably be used. DUSTY is not set up to allow 
dust opacities with radial dependence so only the limiting 
cases with or without ice mantles can be compared. To 
test the sensitivity of the derived envelope parameters on 
the adopted set of dust opacities, the analysis is repeated 



using the bare grain opacities presented by Osscnkopi 



& Hcn aing (1994 ) (column 2 in their Table 1; hereafter 
OH2). From the radial brightness distributions the same 
range of a and AR as when using OH5 is obtained. The 
optical depth, however, is significantly reduced by about 
a factor of two and so is the total mass of the envnlr 



In addition to the analysis of the continuum emission 
the observed molecular line emission is useful in constrain- 
ing the physical properties of the envelope. Traditionally, 
CO and CS line emission have been extensively used for 
this purpose and are adopted here to test the validity 
of the best fit model obtained from the dust anal ysis . 
Using the radiative transfer code presented in Sect. |3.2| , 
the total CO and CS abundances relative to H 2 obtained 
for the best fit model presented in Table f| are ^4xlCU 5 
and ~3xl0 -9 , respectively. In the x 2 -analysis only the 
velocity-integrated intensities in the lines were used. These 
values are within a facto r of about two of what is com- 
monly derived for YSOs ( |van der Tak et al. 2000b| ). The 
abundances in combination with the quality of the fits 
(Fig. [|; see also Sect. ||), in particular the ratios among 
various transitions which are sensitive to the gas temper- 
ature and density, are reassuring and further strengthen 
the adopted physical model. 

In the modelling of the molecular line emission the 
gas temperature is assumed to follow that of the dust. In 
models which self-consistently treat the energy balance the 
gas temperature is generally lower than that of the dust 
i n the outer regions due to imperfect gas-gra i n coupling 



lope 



Adopt ing the model parameters as derived in Tahle Hie et al. 2000a|) . To test the effects of a departure of the gas 



(ICeccarelli et al. 1996| ; |Doty fc Neufeld 1997] ; |Ceccarell 



tioo=4.5, a=1.7, and Ai?=250, but using the OH2 opac- 
ities increases the flux at 60 /im by ~20% compared with 
the OH5 model. The 100 /Ltm flux is not significantly af- 
fected whereas the 1.3 mm flux is roughly twice that of 
the OH5 flux. A model with dust properties varying with 
radius (OH2 in the inner warm parts and OH5 in the outer 
cool parts) would thus serve to improve the fit to the SED. 
Since the bulk of the mass (99.8%) is at low temperatures 
where the grains are coated with ice mantles, the model 
parameters obtained with the OH5 opacities were used in 
the chemical analysis. 



temperature from that of the dust due to gas-grain de- 
coupling in the outer regions, the dust temperature was 
scaled by a constant factor. For Tgas^;0.7xTdust the en- 
velope becomes too cool to fit the observed line intensity 
ratios. Thus, the gas temperature appears to follow that 
of the dust within ^30% in the region probed by the CO 
emission. 

The static envelope model fails to explain the details 
of the individual spectra and the potential of the molec- 
ular emission to constrain the velocity fields will be in- 



vestigated further in Sect. 4.3. The derived abundances 
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Fig. 4. Best fit CO and CS models using a static envelope (dotted line) and a Shu- type collapsing model (full line), 
overlayed on the observed spectra (histograms) (see text for details). The calibration uncertainty in the intensity scale 
for the observed spectra is ~15%. 
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are not very sensitive to the adopted value of a, within 
the limits derived from the dust radiative transfer model. 
Similarly, increasing the outer radius by a factor of two 
only marginally affects the line intensities. The abun- 
dances derived for a wide variety of molecular species are 
further presented in Sect. |5|. 

Thus, the emerging picture from the analysis is that 
the envelope around IRAS 16293-2422 indeed has a re- 
gion of dense and hot gas inside a radius of ^2xl0 15 cm 
(150 AU, 1"), with temperatures decreasing to ~10K at 
-10 17 cm (8000 AU). 

4.3. Infall model 

Current theories of star formation state that protostars 
are formed from the gravitational collapse of cloud cores 
consisting of gas and dust (e.g., Shu et al. 1987). There is 



now growing evidence that low-mass protostars have parts 
of their circumstc llar material in a state of collapse (e.g., 
Myers et al. 200C|) and an obvious extension to the previ- 
ous analysis, which uses a static envelope with a density 
distribution described by a single power-law, is to attempt 
to reproduce the results with an infall model. That the cir- 
cumstellar envelope around IRAS 16293-2422 is in a state 
of collapse has been inferred previously from modelling 
of m il limetre CS observation s ( Walker et al. 1986 ; Zhou 
1995 ; Narayanan et al. 1998 ), although such conclusions 



based on low spatial resolution data have been questioned 
by |Mcntcn et al. (1987 ). Furthermore, Ceccarelli et al. 
(2000a, b) have suggested an infall model based upon a 
physical-chemical model. Here, the observed continuum 
and molecular line emission will be analyzed using the 
well known collapse model presented by |Shu (1977 ). 

In the Shu inside-out collapse model the self-similar so- 
lution is presented in terms of the dimensionless variable 
x—r/at, where r is the radial distance scale, and charac- 
terized by the isothermal speed of sound, a, and the time 
after onset of collapse, t. The location of the collapsing 
wave front at any instant t is described by r c —at. The 
density p and velocity u have the form 



. . a{x) 
PM) = 47rGt2 , u{r,t) = av(x), 



(10) 



where G is the gravitational constant and a(x) and v{x) 
are tabulated by 3hu (1977). In the static part of the en- 
velope {x>l) 



a(x) 



v{x) = 0. 



The asymptotic behavior as x— >0 of the solution is 



a(x) = (S) 1/2 ' v{x) 



2mo 



1/2 



(11) 



(12) 



where mo is equal to 0.975 for this particular solution. 
Given the limited resolution provided by the JCMT at 450 
and 850 microns it is not surprising that the observations 
are well described by a single power-law with a in the 



range ~1. 5-2.0 which are the two extremes obtained from 
the Shu-model. The estimated value of a will depend on 
the location of the collapsing wavefront, r c . In comparison, 
lorgcnsen et al. (2002| ) derived values of a for a large 
number of class protostars and found values typically in 
the range ~1. 3-2.0 . 

The input parameters to DUSTY are the same as for 
the single power-law models (see Table ||) . In addition, the 
envelope size was fixed to a radius 5000 AU. Making the 
envelope larger will produce increasingly worse fits to the 
radial brightness distributions obtained from the SCUBA 
observations. The sensitivity of the two adjustable param- 
eters a and r c =at in the modelling is shown in Fig. || (left 
panel) where the observational constraints used are the 
SED and the SCUBA 450 /xm radial brightness distribu- 
tion. The best fit model is obtained using a~0.6kms _1 
and r c ~2xl0 16 cm, putting the age at ~10 4 yr. 

The mass accretion rate can be estimated from 



M = 



moa 
G 



(13) 



where G is the gravitational constant. For the best fit 
model Af^5xl0 -5 Mq yr _1 is obtained. Such a high ac- 
cretion rate is able to account for the source luminosity. 
The total mass contained in the envelope within 5000 AU 
is ~3.3Mq, consistent with the estimate obtained from 
the static envelope model within the same radius. 

The relative success of the dust modelling using a static 
envelope, with a single power-law to describe the density 
structure, makes it hard to discriminate between the two 
models in the present analysis. In Fig. the density and 
velocity structures obtained from the best fit Shu model 
are presented and, for the density, compared with results 
from the static envelope model. The largest discrepancy 
occurs at small radial distances where the collapsing en- 
velope model predicts about a factor of two to three lower 
densities. We stress the observational data set used in 
the dust modelling is not directly probing this region. At 
larger radii the model is better constrained and the two 
models agree well. It should be noted that the best fit 
single power-law model gives slightly better reduced \ 2 
values for the combined set of observations. However, the 
molecular data provide further constraints since they have 
the potential to probe the large scale velocity field. 

In Fig. [I], spectra of CO and CS line emission as ob- 
served with the JCMT are presented. Lines which are opti- 
cally thick, like CS, 13 CO and C 18 (J=3->2), show a dis- 
tinct, narrow, absorption feature near the stellar velocity. 
This feature is due to effective self- absorption in the outer 
cool parts of the envelope. Also, the degree of the asym- 
metry in the line profiles increases with the optical depth 
in the lines. In the optically thin lines the self-absorption 
feature disappears and the lines are well described by a 
single gaussian profile. The width of the self-absorption 
feature constrains the turbulent velocity to ^0.3 km s -1 
in the outer envelope. For simplicity this value is adopted 
throughout the envelope. A turbulent velocity component 
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Fig. 5. Properties of the circumstellar envelope around IRAS 16293-2422 obtained from the dust modelling. Shown 
are the dust temperature (left) and density (middle) structures for best fit model (full line) assuming a single power- 
law distribution of the density. In the temperature panel, an optically thin prediction using an opacity law k v qoj^ 
with (3=1 (dotted line; Tdocr -0 - 4 ) and (3=2 (dashed line; Tdtxr -0 33 ), is also shown. In the density panel the best fit 
Shu- type collapsing envelope model is also shown (dash-dotted line). The velocity structure obtained from the Shu 
infall model is shown in the right panel. 
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Fig. 6. x 2 " m aps showing the sensitivity of the Shu-collapse model to the adjustable parameters in the modeling of 
the dust continuum emission (left) and the CS (CS and C 34 S transitions; centre) and CO ( 13 CO, C 17 0, and C 18 
transitions; right) line emission compared to observations. Contours are drawn at x min +(2.3, 4.6, 6.2) indicating the 
68% ('lcr'), 90%, and 95% ('2<t') confidence levels, respectively. The number of observational constraints used, N, are 
also shown. The quality of the best fit model can be estimated from the reduced chi-squared statistic Xred = X m in/(-^~ 2). 



varying with radius is beyond the scope of this article, see 
however Stark et al. (2002, in prep.). 

The observed CS emission (including that from C 34 S) 
is analyzed using the abundances derived from the static 
envelope model presented in Sect. (see also Sect. ||). 
From only the integrated intensities it is possible to con- 
strain both a and r c as is shown in Fig. a (middle 
panel). The model spectra are presented in Fig. |3 together 
with the observations. The fit to the integrated intensi- 
ties is worse than that obtai ned from the single power- 
law model presented in Sect. 4.2, however. In particular 
the C 34 S(J=7— >6) line is poorly reproduced in the infall 
model. Using also the line profiles as constraints we find 
that a model where the collapsing wavefront is located 
at — 1.0xl0 17 cm and the value of a is ^0.9kms _1 best 
reproduces the observations. 

Analyzing the CO emission ( 13 CO, C 18 0, and C 17 0) 
gives yet another set of estimates. At first the abun- 



dances derived from the static envelope models presented 
in Sect. [T| (see also Sect. |^) are adopted. From the anal- 
ysis of the integrated intensities shown in Fig. 6j (right 
panel) and the line profiles presented in Fig. | a Shu- 
model with a~1.0kms _1 and r c ~1.2x 10 17 cm reproduces 
the CO observations well, in excellent agreement with the 
CS modelling. The derived dynamical age of the system 
is — 3xl0 4 yr. In contrast with the CS modelling, the fit 
to the integrated CO intensities is equally good as ob- 
tained for the single power-law model. As discussed in 
Sect. ||, the CO abundance obtained from the static en- 
velope model is about a factor of 2—3 lower than what is 
typically observed for interstellar gas. If instead the CO 
abundance relative to H2 is assumed to be the 'standard' 
interstellar value of lxlO~ 4 and the standard isotopic 
ratios are assumed ( [CO/ 13 CO]=60, [CO/C 17 O]=2500, 
and [C 18 0/C 17 0]=3.9) the estimate of a is -0.75 and 



F. L. Schoier et al.: Does IRAS 16293-2422 have a hot core? 



13 



r c ^3x 10 17 cm. However, the quality of the fit becomes 
worse in this case. 

The results obtained from the best fit Shu-model are 
presented in Fig. [| overlayed onto the observed CO and 
CS spectra. The integrated intensities and, to some ex- 
tent, line profiles can be modelled with the spherically 
symmetric infall solution. The details of the spectra will, 
however, intricately depend on the adopted geometry, ve- 
locity fields, and chemical gradients. The velocity field, 
in particular position-velocity maps of the source, form a 
stringent test of the dynamical models. For IRAS 16293- 
2422 it will likely be necessary to include a rotational com- 



poncnt to the velocity f ield (Mcntcnct al. 1987; Zhou 1995; 
Narayanan et al. 1998) . 



Other estimates of the infall radius, based on analy- 
sis of molecular line emission, range between about 5 and 
15xl0 16 cm (|Walker et al. 1986|; ^hou 1995|; |Narayanan 



et al. 199S; |Ceccarelli et al. 2000a[) in excellent agree- 
ment with the values obtained here from analysis of CO 
and CS emission. However, there appears to be a discrep- 
ancy between the dust and molecular line analysis, possi- 
bly reflecting the fact that the simple Shu-collapsing core 
model is not fully adequate to describe the state of the 
infalling material and/or that some of the CO and CS 
emission is associated with the outflow and surrounding 
cloud. Moreover, gradients in the molecular abundances 
will affect the parameters derived. 

5. Molecular abundances 

The basic envelope parameters derived from the dust ra- 
diative transfer modelling performed in Sect. ^, in par- 
ticular the density and temperature distributions, are 
used as input for the Monte Carlo modelling of the 
molecular line emission. The static power-law model is 
adopted; the abundances obtained with the best fitting 
infall model generally differ by no more than ^25% for 
a constant abundance model. However, models where 
a drastic enhancement in the abundance is introduced 
('jump-models') require 2—3 times larger abundances in 
the inner hot part for a Shu-type collapsing model, reflect- 
ing the significantly lower density compared to the static 
power-law model in this region (Fig. |^). Changing the en- 
velope parameters describing the static power-law model 
within the accepted range of values (Fig. ^ typically af- 
fects the abundances obtained for the best fit model by 
less than ±25% for constant abundance models. In 'jump- 
models' the effect on lines which are sensitive to the condi- 
tions in the innermost dense and hot regions can be higher, 
up to ±50%. 

5.1. Constant abundance models 

The abundances are initially assumed to be constant 
throughout the envelope. For simplicity the observed lines 
are assumed to be broadened (in addition to thermal line 
broadening) by microturbulent motions only. The micro- 
turbulent velocity is set equal to 2kms _1 throughout the 



envelope typically producing lines ~4kms _1 (FWHM) 
wide (see Sect. ||). The presence of a global velocity field, 
e.g., infall, outflow, or rotation, would serve to reduce the 
optical depths of the line emission, so that the abundances 
presented in Table [sj are strictly lower limits. However, 
such effects will not significantly increase the inferred 
abundances since they are largely derived from optically 
thin lilies. The observational constraints used in the mod- 
elling are the total velocity-integrated line intensities. 

In many cases the observations are well reproduced as- 
suming a constant abundance throughout the envelope, as 
seen from their reduced x 2 ~l- I n the cases where the fits 
are good the derived abundances are generally consistent 
with typical values found in quiescent molecular clouds. In 
addition, the isotopic ratios of 18 0/ 17 0~3.9 determined 
from CO observations and 32 S/ 34 S~25 from CS observa- 



tions, agree well with interstellar values (Wilson fc Rood 



1994 ). A notable exception is the relatively low abun- 
dance, ~7xl0 -11 , derived for HNC. The abundances de- 



rived by Blake et al. (1994 ) and van Dishoeck et al. (1995 ) 
agree surprisingly well with the new, more accurate, esti- 
mates presented here (Table |]), typically within a factor of 
~2. Those abundances were derived from statistical equi- 
librium equations assuming a constant temperature and 
density. The agreement indicates that the adopted values 
were representative of the region from which most of the 
submillimetre emission arises. 

For the main isotopes of HCN and HCO + the emis- 
sion is highly optically thick and the models are rela- 
tively insensitive to the molecular abundance. The abun- 
dances of these molecules were instead estimated from the 
rarer isotopomers assuming a standard isotope ratio, i.e., 
12 C/ 13 C=60. The abundances obtained in this way fail 
to account for all of the observed flux in the HCN and 
HCO + lines, as evidenced by their high reduced x 2 -values 
in Table H Adopting the best fit Shu-infall model derived 
from CO and CS observations and presented in Sect. 4.2 
introduces a large-scale velocity field which reduces the 
line optical depths and increases the line intensities thus 
improving the fit to observations. Material in the outflow 
can also contribute to these lines. 

The observed transitions are only in LTE throughout 
the envelope for abundant molecules like CO and OCS, in- 
cluding their isotopomers observed here. For less abundant 
species, where collisional excitation is less efficient, depar- 
tures from LTE are found. For example, the level popula- 
tions of common molecules like CS and H2CO are in LTE 
out to ~1— 2xl0 16 cm. For most molecules, populations of 
the observed lines are in LTE within ^2xl0 15 cm. 

From Table | it is also evident that the line emission 
from several molecular species is not fitted well, in par- 
ticular molecules where the emission probes a large radial 
range, e.g., in the case of H2CO and CH3OH. In addition, 
the isotopic ratios derived in many of these cases are far 
from their interstellar values and what is commonly de- 
rived for these kind of objects. Typically, the model inten- 
sities from lines sampling the inner parts of the envelope 
are too low compared to observed values, whereas the op- 



14 F. L. Schoier et al.: Does IRAS 16293-2422 have a hot core? 

Table 5. Derived abundances using a constant molecular abundance fx relative to H2 throughout the envelope. 





Molecule 




/ 


X 




2 a 
Xrcd 


N h 


rold c 

Jx 




P d 

[K] 


rp c 
-C-'max 

[K] 


13 HO 

V V ./ 


Ldl UU11 1I1LJ11U.A.1U.C 


U . J 




1 0" 


7 




■i 




1.6 


X 


10" 


T 


1 5 Q 


31 7 


c 18 o 




6.2 




10" 


8 


0.5 


2 


1.0 


X 


10" 


7 


15.8 


31.6 


r 17 n 

V KJ 




1 

1.0 


X 


1 n~ 


8 


n 1 

U.l 



z 


3.8 


X 


10" 


8 


1 9 


39 A 


lluU 


1U1 illy 1 lUil 


1.4 




1 0" 


1 


1 n q 


Q 
O 


1.8 


X 


10" 


T 


95 7 


42.8 


1 1 v V 7 




9 4 


X 




11 


n 1 

u. 1 


O 


7.5 


X 


10" 


12 


9^ n 


41 fi 


Hr 18 n+ 






X 


1 nr 


12 




1 

1 


3.5 


X 


10" 


12 






1 J V. \J 




1 3 




1U 


11 


u . 


2 


1.5 


X 


10" 


11 


90 7 

ZO. 1 


51 Q 




c y diiu^tjii 


O.U 


X 


1 n~ 


11 


1 9 


A 


1.0 


X 


10" 


10 


1U.O 


7 

oz. / 


HCN 


livrlrncpn pvn n irlp 

LL y CfciC±± y CL±11V_1.C 


1.1 




10" 


1 


2.3 


3 


1.9 


X 


10- 


9 t 

1 


25.5 


42.5 


H 13 CN 




1.8 




10" 


11 


6.2 


2 


1.4 


X 


Hr 


11 


24.9 


41.4 


HC 15 N 




7 fi 




1 n~ 

J-U 


12 






7.0 


X 


10" 


12 


24.8 


24.8 


DCN 




1 

D . O 




1 n~ 

J-U 


11 


. 


2 


2.5 


X 


10" 


11 


9D Q 
zu. y 


52.1 


HNC 


11 V Lll Utcil loUcy clillUC 


Ci Q 

V) . i7 




1U 


11 


n 9 


Q 
O 


1.5 


X 


10" 


10 


26.1 


43 5 


HN 13 C 




5.2 


x 


10" 


12 






2.5 


X 


10" 


12 


25.1 


25.1 


DNC 




4.2 




10" 


12 


3.0 


2 


5.0 


X 


10" 


12 


22.0 


36.6 


HC-)N 


pvfl n oa ppt vl pti p 


1.5 




10" 


10 


1.8 


3 


2.5 


X 


10" 


11 


131.0 


177.3 


V_^ll3 ' 


tyic»t n^rl *"M 7 c* tii /"i 
lilc Lll V icy dlllUc 


1.4 




1 n~ 

J-U 


10 




7 


1.5 


X 


10" 


1(1 




957 Q 
zj 1 . y 


m\ v.. v 7 


IbUcyeLillc dclLl 


1 


X 


1 nr 

J-U 


10 


1 n 




1.7 


X 


10" 


10 


1 1 9 fi 


1 9fi fi 




c Lrlly lly 1 


2.1 


y 


1 n~ 

1U 


10 


u . 


4 


2.5 


X 


10" 


10 


25.1 


fi9 Q 
uz. y 










1U 


11 


n n 


2 


4.5 


X 


10" 


11 


9D 7 

ZU. 1 


90 7 

ZU. / 


C3H2 


cy ClULliUJJciiy HVlcilc 


1 fi 


X 


1 n~ 


11 


A 1 


ft 
u 


3.5 


X 


10" 


11 


1 Q f? 
iy.j 


8fi Q 
ou.y 




m t \\~%.t\ a ^dti/I on 

nic oily 1 d/Oc Ly iciic 


1.4 


y^ 


1 0~ 

J-U 


9 




U 


6.5 


X 


10" 


10 


74 fi 


1 fi3 9 


u-ii2v^'v-' 


TAVm !S 1 f \ OrH7n 
lUi illdlLlciiy Uc 


U ■ 3 




1 n~ 

J-U 


10 


4.1 


q 


5.2 


X 


10" 


10 


21.9 


1 74 n 


n H»rn 

p-n2V^w 




r i 1 
•J ■ ± 


X 


1 nr 

J-U 


1(1 


9 Q 





1.8 


X 


10" 


10 


91 n 

Zl.U 


94PI 7 


_H 2 13 CO 




1.2 




10" 


11 


3.7 


4 










21.7 


61.3 


p_H 2 13 CO 




0.5 




10" 


12 


3.6 


2 










51.1 


98.4 


n-HDCO 




1.0 




10" 


10 


0.1 


2 


7.3 


X 


10" 


11 


52.4 


56.3 


1 I 1 1 ^ V V ) 




■ ± 


X 


1 n~ 

J-U 


11 


. u 


9 
z 


2.5 


X 


10" 


11 


ou.o 


fi9 7 
uz. / 




TTl 0"f~ Vl ' 1 11 

111c LlldllUl 


1 7 


y^ 


J-U 


9 


0. 1 


9^ 
z.> 


4.4 


X 


10" 


9 


1 5 5 
10.0 


1 87 fi 






0.0 


X 


1 nr 

J-U 


10 


9 


9 
z 


1.8 


X 


10" 


10 


00. u 


1 fin n 


HrnnH 

1 1 V. V y V y I 1 


tcwyywc of in 
lUliillC elClLl 






1 n~ 

1U 


10 






<3.0 


X 


10" 


10 






v_^ii3 1 1 w 


dec LcLlVlciiy Vic 




y 


1 n~ 

1U 


11 






<1.0 


X 


10" 


10 








H l yyi ot n\r 1 n cv 
Lllillc Lrliy 1 c Lilc! 




y 


1 n~ 

1U 


10 






<2.0 


X 


10" 


9 






1 1 V. V 7 V 7 v../ 113 


vri t n^r 1 ti"iT"Tti uro 
lilc Lily 1 1U1 illelLc 






1 0~ 

1U 


9 




















cell UUil IIILJIIUO Lllllclc 






1 0~ 

J-U 


9 




3 


1.1 


X 


10" 


9 


35 3 
00.0 


fi5 8 


c 34 s 




1.2 




1 n~ 

J-U 


1(1 


n q 


3 


5.0 


X 


10" 


11 


34 7 


fi4 8 


so 


qii lfiir mnnn^riolp 

O U.l± Lll. ± 1 J_ \_f L 1 UA 1 LI c 


4.4 




10" 


9 


2.4 


() 


3.9 


X 


10" 


9 


35.0 


87.5 




ciiItiit" Hi rw ino 
Lill HI UlL)Jill_lc 


f-"i 9 

V) . z 


y 


1 n~ 

1U 


10 


fi 1 

U. J- 


1 n 


1.5 


X 


10" 


9 


19.2 


89 5 


or 1 


caxbonyl sulfide 


7 ("1 


X 


1 n~ 


9 


u . 


9 
z 


7.1 


X 


10" 


9 


1 99 ft 


9*^7 n 
Zo / .u 


1 IV O 


thioformyl ion 


9 1 


X 


1 nr 


11 


t .u 


9 
Z 


2.0 


X 


10" 


11 


7 


7^J 7 


o 13 cs 




6.7 


X 


10" 


10 


5.6 


2 


2.8 


X 


10" 


10 


99.5 


122.2 


OC 34 S 




1.6 


X 


10" 


9 


1.8 


3 


3.2 


X 


10" 


10 


99.8 


147.7 


H 2 S 


hydrogen sulfide 


1.6 


X 


10" 


9 




1 


1.5 


X 


10" 


9 


84.0 


84.0 


0-H2CS 


thioformaldehyde 


2.0 


X 


10" 


10 


8.8 


6 


1.1 


X 


10" 


10 


58.6 


209.1 


p-H 2 CS 




1.8 


X 


10" 


10 


2.9 


4 


5.7 


X 


10" 


11 


46.1 


98.8 


SiO 


silicon monoxide 


5.9 


X 


10" 


11 


2.2 


4 


1.0 


X 


10" 


10 


31.3 


75.0 


29 SiO 




8.8 


X 


10" 


12 


5.5 


3 


5.0 


X 


10" 


12 


31.0 


74.3 



a The reduced x 2 of the best fit model. Based on the x 2 - ana lysis the derived abundances are accurate to about 20—30% 
when the fit is good, i.e., x 2 c d~l- 

b The number iV of inde pendent observatio nal c onstraints used in the mod elling. 
c Abundance derived by 



Blake et al. (1994) and van Dishoeck et al. (1995) using a simple excitation analysis. 



Lowest energy of the upper level involved in the transitions used as constraints in the modelling. 
c Highest energy of the upper level involved in the transitions used as constraints in the modelling, 
f Derived from optically thin isotope(s) assuming standard isotopic ratio(s). 



F. L. Schoier et al.: Does IRAS 16293-2422 have a hot core? 



15 



posite is true for the lines probing the outer part of the Table 6. 
envelope. An obvious explanation is that a steep gradient a jump i 
is present in the abundances of these molecules. T=90K. 



Various molecular abundances derived using 
i their fractional abundance, introduced at 



5.2. Jump models 

In order to improve the quality of the fits, models with a 
jump in the molecular abundances are considered. A jump 
is introduced at the radius in the envelope where the tem- 
perature reaches 90 K, at which point the ice starts to 
evaporate from the grain mantles. In our models, this oc- 
curs at ^2xl0 15 cm or 150 AU. The free parameters in 
the modelling are then the fractional abundances in the 
inner (/i n ; T>90K) and outer parts of the envelope (/ ou t; 
T<90K). For most species, the isotopic abundance ra- 
tios are assumed to be fixed to the standard interstellar 
values to increase the number of constraints used in the 
modelling. The results from the excitation analysis are 
presented in Fig. [?] for several species. The reduced \ 2 
for the best fit models are generally very good, ~1, and 
significantly better than the constant abundance models. 
Typically, a jump of ^100 in abundance is derived, leading 
to abundances that are significantly higher than found in 
quiescent molecular clouds and comparable to those found 
in the prototypical hot core in Orion (Sect. ||). 

The jump models can only be applied to species for 
which a significant number of lines are observed covering a 
wide range of excitation conditions. These include H 2 CO, 
CH 3 OH, CH 3 CN, H 2 CS, SO and S0 2 (see the columns of 
E max and E min included in Table |). For HC 3 N and OCS, 
only lines from highly-excited levels have been observed, 
so that for these molecules the values of / ou t in the outer 
envelope are poorly constrained. For most simple linear 
rotors such as HCN, HCO + , CN, however, the observed 
lines arise from levels below 90 K, so that no informa- 
tion on the inner warm part is obtained. The only excep- 
tion is SiO, where the combination of many 28 SiO and 
29 SiO lines allows a jump to be inferred. Molecules such 
as HNCO, CH2CO and H2S for which only one or two 
lines are observed and where the emission mainly probes 
hot gas (Table |^) only the inner part of the envelope were 
modelled. The significantly higher abundances obtained 
(Sect. H) compared with the constant abundance mod- 
els illustrates the point that orders of magnitude higher 
abundances can be derived if the emission is assumed to 
originate only from the inner warm region. 

The observations analyzed in this paper do not in- 
clude the lowest rotational transitions probing the coldest 
outer parts; thus, so-called 'anti-jump' models, in which 
the abundances are decreased below a certain temperature 
due to freeze-out, cannot be tested, except for the case of 
CO. There are also some molecules, e.g., C3H2, for which 
the jump-models give a worse x 2_n t than the constant 
abundance models; such molecules are good candidates 
for the 'anti-jump' models if lower transitions are avail- 
able. In the following, a few individual cases are described 
in more detail, before discussing the general results. 
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a Abundance in the inner, dense, and hot part of the envelope. 
b Abundance in the cooler, less dense, outer part of the 
envelope. 



5.3. CO 

Carbon monoxide, CO, is difficult to destroy but relatively 
easy to excite through collisions, even in the outer low- 
density and cold part of the envelope. Brightness maps of 
the 12 CO molecular line emission associated with IRAS 



1629 3-2422 reveal a complex structure ([Walker et 
1988 ) indicating 



al 



two 



bipolar outflows. Intcrfcrometric 
BIMA observations suggest that the 13 CO(J=1^0) emis- 
sion is also associated with the outflow to some extent 
(Schoier et al. 2002a, in prep.). Single-dish observations 
of higher transitions of 13 CO show no direct evidence for 
tracing the outflow based upon the shape of their line 
profiles (see Fig. |]). This is also the case for lines from 
the less abundant C 17 and C 18 molecules. Thus these 
lines can possibly be used to trace the CO content in the 
circumstellar envelope. 

Assuming a constant abundance throughout the enve- 
lope the radiative transfer calculations give abundances 
of 6.5xl0~ 7 , 6.2xl0~ 8 and 1.6xl0~ 8 for 13 CO, C 18 0, 
and C 17 0, respectively. The derived C 18 0/C 17 ratio 
is 3.9, in good agreement with typical interstellar val- 
ues QWilson fc Rood 1994] ). [l0rgensen et al. (2002| ) 
derived i8 0/ 1Y 0=3.6±l.l for their survey of proto- 
stcllar objects. The 12 CO abundance is estimated to 
be 4.0xl0~ 5 from C 17 assuming the terrestrial ratio 
12 CO/C 17 O=2500. This is in excellent agreement with the 
value of 3.9xl0 -5 estimated from 13 CO using the inter- 
stellar value 12 CO/ 13 CO=60. In addition, the upper limits 
obtained for the 13 C 17 0( J=2->1) and 13 C 18 0( J=3^2) 
line emission are also consistent with these values. In all, 
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Fig. 7. x 2 -maps of various molecular species introducing a step function in abundance at T gas =90K. Contours are 
drawn at Xmin+(2-3, 4.6, 6.2) indicating the 68%, 90%, and 95% confidence levels, respectively. The number of ob- 
servational constraints used, N, are also shown. The quality of the best fit model can be estimated from the reduced 
chi-squared X r 2 cd =X^ in / '(JV "-2) . Fixed isotope ratios of 28 SiO/ 29 SiO=20, OCS/OC 34 S=20 and OCS/O 13 CS=60 were 

assiimprl 



F. L. Schoier et al.: Does IRAS 16293-2422 have a hot core? 



17 



the consistency of the derived values and quality of the 
model fits (Fig. ^ and Table ||) are reassuring. Due to the 
high optical depths in the observed 13 CO lines, their pro- 
files tend to be somewhat broader than those from the less 
abundant isotopomers, and their intensities are less sensi- 
tive to the assumed abundance. The observed lines are all 
in LTE, indicating that the derived abundances are not 
sensitive to the adopted set of collisional rates. 

The total inferred CO abundance is about a factor of 
two to f our lower than the va lue of 8 x 10~ 5 found in dark 
clouds ( Frerking ct al. 1982 ) and 2 x 1CP 4 in warm re- 



gions ( Lacy et al. 1994 ). A plausible explanation is that 
CO freezes out in the cool external parts of the envelope. 
To simulate this situation, a 'jump' model with an abrupt 
decrease in the CO abundance at 20 K was introduced; 
this is a characteristic temperature below which pure-CO 
ice can exist ( Bandford fc Allamandola 1993| ). To com- 
pensate for this freeze-out, the CO abundance in regions 
above 20 K needs to be raised. In Fig. || the result of vary- 
ing the CO abundance in the inner and outer parts of the 
envelope is presented. A maximum allowed value for the 
CO abundance of gas above 20 K is —7xl0~ 5 while at 
the same time the abundance in the outer cooler envelope 
needs to be lowered to <2x 10 -5 , i.e., a depletion of about 
a factor four. However, in the present analysis a constant 
abundance model is equally probable. The molecular line 



modelling performed in Sect. 4.2 suggests that in the en- 
velope around IRAS 16293-2422, T gas >0.7xT dust . This in 
turn means that it is not possible to explain the appar- 
ently low CO abundance with a significant decoupling of 
the gas temperature from that of the dust. 

For the similar modelling of the larger sample of class 
and I objects, |j0rgensen et al. (2002| ) found that the CO 
abundance in general was lower for the class objects than 
the class I objects (average CO abundances of respectively 
2.0 xlO" 5 and 1.2xl0" 4 ). Further it was found that abun- 
dance jumps at 20 K of more than a factor 3 could be ruled 
out in most cases and that constant fractional abundances 
over the temperature range covered by the CO rotational 
lines provided good fits. This lead to the suggestion that 
CO in the class objects could be trapped in a porous 
ice matrix with H2O from which it does not fully evapo- 
rate until at temperatures of ~ 60 K. More observational 
constraints on isotopic CO ( J=l— >0) as well as higher— J 
CO lines are needed to verify if CO is frozen out onto dust 
grains only at the lowest temperatures or if a substantial 
fraction of CO evaporates more gradually up to —90 K, as 
suggested by recent experiments of CO-H2O ice mixtures 
QCollings et al. 2002[ ). 



In Fig. |] the line profiles obtained from the constant 
abundance model of the CO emission are presented. The 
asymmetry present in some of the observed spectra is not 
possible to model using a static envelope and requires the 




2x10" 



Fig. 8. CO jump-model assuming an abrupt change in 
abundance at 20 K using the integrated intensities of the 
observed 13 CO, C 18 0, and C 17 line emission as con- 
straints. Contours are drawn at x, 2 nin +(2.3, 4.6, 6.2) in- 
dicating the 68%, 90%, and 95% confidence levels, respec- 
tively. The number of observational constraints used, AT, 
are also shown. The quality of the best fit model can be es- 
timated from the reduced chi-squared Xred = Xmin 

/(N-2). 



5.4. H 2 CO and CH 3 OH 

H 2 CO and CH 3 OH are two particularly useful molecules 
to model in detail, because of the relative complexity 
of the energy level structure of asymmetric rotors and 
symmetric top molecules. For example, for H2CO dipole 
transitions between various -£T p -ladders are not allowed 
while collisional transitions are, making these transitions 
good probes of the temperature structure. In addition, 
the wealth of available lines in the submillimetre regime 
makes them ideal to probe various parts of the envelope 



(Mangum & Wootten 


1993 


Ceccarelli et al. 2000b 


). It is 



van Dishoeck et al. 1993 ; 
:lear from Tabic |6j that the 
reduced \ 2 values are considerably lowered in the jump- 
model for an abundance contrast of —60— 100 for H2CO 
and -50 for CH 3 OH. The H 2 CO and CH 3 OH abundances 
in the cool outer part of the envelope are in line with values 
in quiescent molecular clouds, whereas those in the inner 
envelope are of the same order of magnitude as derived 
for some hot cores (Sect. |). Typical CH 3 OH and H 2 CO 
abundances in interstellar ices are <2xl0 -6 and — 10~ 6 



(Dartois et al. 1999; Keane et al. 2001), illustrating that 
the drastic increases could be explained by ice evaporation 
from dust grains above 90 K. For high-mass protostars, 
van der Tak et al. (2000a) found a similar CH 3 OH abun- 



presence of a global velocity field, as discussed in Sect. 4.3 



dance jump for some sources, but no evidence for signifi- 
cant jumps in the H2CO abundance. The H2CO/CH3OH 
ice abundance ratio is sensitive to the atomic hydrogen 
density in the ice-forming region, which could be different 
for low- and high-mass protostars. 
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When a significant number of transitions with vary- 
ing excitation conditions are observed the data have the 
potential to determine the characteristic temperature at 
which the majority of the molecules are evaporated. The 
sensitivity of the data presented here to the adopted jump- 
temperature can thus be tested for H2CO and CH3OH. 
The H2CO models are not very sensitive to the jump 
temperature until it drops below about 40 K at which 
point the reduced x 2 increases fast. It is found that a 
temperature of ^50 K gives the best fit indicating that 
formaldehyde starts to evaporate at temperatures be- 
low 90 K (see also Ceccarelli et al. 2001). In the case of 
methanol (CH3OH), the best fit is obtained for ^90 K 
whereas temperatures below ~50K give poor fits, indi- 
cating that this molecule evaporates mainly at ^90 K. 

The ortho-to-para ratio of H2CO is not well con- 
strained. Based on the x 2 -analysis, its lower limit is ~0.9, 
with a best-fit value of 2.5. A further complication is 
that the ortho-to-para ratio may vary, e.g. according to 
temperature, through the envelope. In the present anal- 
ysis it is not possible to confirm this. Fixing the ortho- 
to-para ratio to 2.5 throughout the envelope increases 
the number of constraints used in the modelling and en- 
ables 'jump-models' for H2 13 CO and HDCO to be in- 
ferred. From the isotopic ratio, the derived 12 C/ 13 C-ratio 
is ~100 with considerable uncertainty and fully consis- 
tent with the interstellar value of 60 adopted elsewhere 
in this paper. From deuterated formaldehyde (HDCO) it 
is possible to estimate the D/H ratio and a value ~0.3 
is derived. This is significantly higher than what is ob- 
tained from DCO+/HCO+, DCN/HCN and DNC/HNC 
(see Table ||), and suggests a different scenario for H2CO. 
The high degree of deuterium fractionation of H2CO, fur- 
ther str engthened by the recent detection of D2C O in this 
source ( Loinard et al. 2000 ; Ceccarelli et al. 2001 ) , is about 
five times larger than values obtaine d towards other low- 
mass protostars ( Roberts et al. 2002 ) and dark clouds like 
TMC-1 and L134N ( |Turner 200l|) but consis tent with es- 
timates for the Orion hot core ( Turner 1990 ). Due to the 
limited number of lines any radial variations in either the 
12 C/ 13 C-rat io or the deuterium fr actionation cannot be 
established. Ceccarelli et al. (2001 ) argue, based on spa- 
tially resolved emission, that H2CO and its deuterated 
counterparts are formed mainly from grain-surface reac- 
tions in a previous cold, dark cloud phase. For HCO + , 
HCN and HNC the degree of deuteration can be explained 



by gas-phase reactions at low temperatures (Roberts et al 
20020 



5.5. Vibrationally excited emission? 

For linear rotors such as CS and HCN, the observed lines 
cover only a limited range in excitation energy. Much 
higher frequency data are needed to probe the abundances 
of these molecules in the warm gas through pure rotational 
lines in the vibrational ground state (e.g., Boonman et al. 
2001). Some lines of vibrationally-excited molecules have 



been detected toward IRAS 16293-2422, however, in par- 
ticular the CS(i>=l; J=7^6) line. Since this line originates 
from a level at about 1900 K above ground, it cannot be 
excited by collisions only. Instead, radiative excitation by 
dust emission must play a role. 

The addition of radiative excitation by dust is straight- 
forward in the Monte Carlo scheme if scattering is as- 
sumed to be negligible, which is the case at the wave- 
lengths of importance here: the CS fundamental vibra- 
tional transition occurs around 8 ^tm. In the absence of 
scattering events only emission and absorption by the dust 
particles need to be considered. The dust is assumed to 
locally emit thermal radiation described by the dust tem- 
perature T(r), according to Kirchhoff's law. The model 
photons emitted by the dust are released together with 
the other model photons and the additional opacity pro- 
vided by the dust is added to the line optical depth. 

The CS abundance of2.5xl0~ 9 derived previously fails 
to account for the observed CS(t;=l; J=7^6) line emission 
by many orders of magnitudes. Even a 'jump-model' with 
the abundance increased by a factor of 100 fails to account 
for all the observed emission. 



Recently, Highberger et al. (200C ) failed to detect any 
vibrationally excited CS emission in lower J-transitions 
towards IRAS 16293-2422 down to T rms ~10mK. One 
possibility is that the CS(t>= l; J=7^6) line was mis- 
identified by Blake et al. (1994 ). Their observations were 
obtained in dual sideband and the CS line was supposed 
to reside in the lower sideband. However, in the up- 
per sideband two H 2 CS lines {Jk~ ,K+ = ^5,6~*^s,s an d 
J K _ jK+ =105,5^-95,4) at 343.202331 GHz coincide with the 
position of the CS line in the lower sideband. To check 
the possibility that these transitions contribute signifi- 
cantly to the observed intensity of the line, a LTE jump- 
model was run using the parameters derived previously 
for H2CS. It is indeed found that these high lying energy 
levels (-E U =418K) are sufficiently excited to account for 
all of the flux observed in this line. 



6. Discussion 

6.1. Does IRAS 16293-2422 have a hot core? 

Hot cores have been defined for high-mass sources as small 
(<0.1pc ), dense (n H , >10 7 cm' 3 ) and wa rm (T>100K) 
regions (Walmsley 1992; Kurtz et al. 200C). The chemical 
signatures of hot cores are high abundances of a wide vari- 
ety of complex organic molecules and fully hydrogenated 
molecules such as H2O, NH3 and H2S. These chemical 
characteristics are thought to arise from thermal evapo- 
ration of ice mantles from the grains close to the proto- 
star, followed by rapid high-temperature gas-phase reac- 
tions for a period of 10 4 — 10 5 yr (Charnley et al. 1992). 



In this scenario, the fully hydrogenated species are 'first 
generation' molecules produced by surface chemistry on 
the grains, whereas complex organic molecules like HC3N, 
CH3OCH3 and HCOOCH3 are 'second generation' pro- 
duced by gas-phase chemistry between evaporated species. 



F. L. Schoier et al.: Does IRAS 16293-2422 have a hot core? 
Table 7. Summary of derived abundances in the envelope of IRAS 16293-2422. 
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a Abundance in warm and dense inner part of the envelope <150 AU in radius around IRAS 16293-2422 

Abun dance in cooler, les s dense outer part of the env elope around IRAS 1 6293-2422 
c From lOhishi et al. (1992), updated using values from Dickens et al. (2000) for L134N pos. C 
van^J2ish£C£k^£l31ake^998 ) (and references cited); otherwise from Sutton et al. (1995) 
Bachiller fc Perez Gutierrez (1997 ) at position B2 assuming CO/H2=10 _ 
f Based on Ehrenfreund & Charnley (2000) assuming H2O(ice)/H2=5xl0 -5 . The ranges reflect the values found for different 

assuming this ice evaporates as HNCO 

-5 



Froml 
From 



sources. The value for HNCO is derived from OCN 



s From Bockelee-Morvan et al. (2000) for comet Hale-Bopp at 1 AU, assuming H2O/H2=5xl0 
Assuming that the observed emission or upper limits apply to the inner warm region only 



To what extent do these physical and chemical character- 
istics also apply to the low-mass object IRAS 16293-2422? 

The dust radiative transfer modelling performed in 
Sect. [| clearly indicates that there is evidence for the pres- 
ence of hot and dense material within ^150 AU of the 



protosiiar. Moreover, this material is most probably m a 



state of collapse towards the protostar. The need for m- 
troducing drastic jumps in the abundances of some 'first 
generation' species like H 2 CO and CH 3 OH at the evap- 



comet Hale-Bopp, which may be representative of inter- 
stellar ices. Although the IRAS 16293-2422 abundances 
are still up to an order of magnitude lower than those 
found in typical ices, they are comparable to those ob- 
served in high-mass hot cores. The high degree of deu- 
terium fractionation measured for these molecules (van 
Dishoeck et al. 1995 ; |Ceccarelli et al. 2001 ) is consistent 



oratior j~tcmpcraturc of the ices (^90 K) further suggests 
that the ice mantles are liberated 111 the hot inner regions. 
Table | compares the abundances f m of these molecules 
with those found in interstellar ices and in high-mass hot 
cores. The table also includes the abundances found in 



with this interpretation, and probably results from a com- 
bination of gas-phase deuterium fractionation in the cold 
pre-collapsing cloud and grain-surface reactions (Ticlent 
1983| ; [Roberts fc Millar 2000| ). 



The high abundances of the sulfur-bearing species of 
~10 -7 may also fit with this scenario. Molecules like SO, 
SO2, OCS and H 2 CS are all predicted to be drastically 
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enhanced by the injection of H2S into the hot core region 
(Charnley 1997), but may also be present in the grain 
mantles themselves. H 2 S is detected toward IRAS 16293- 
2422, but because only a single line is observed, the jump 
model cannot be uniquely constrained. If all H 2 S emis- 
sion is assumed to come from the hot inner region, its 
inferred abundance is ~lxl0~ 7 . This relatively high H2S 
abundance is consistent with values found in a survey of 



massive hot cores (Hatchcll et al. 1998) but an order of 



magnitude lower than the estimates for Orion (Minh et al 
199C0 

The physical characteristics in IRAS 16293-2422 are 
thus consistent with those of a typical 'hot core' as ob- 
served towards many high mass protostars, except for a 
large change in physical scale. The remaining question is 
whether the evaporated ices have driven a similarly com- 
plex organic chemistry in this low-mass protostar. Among 
potential 'second generation' products, HC3N and CH3CN 
are observed in fairly large amounts with some evidence 
for jumps in their abundances. However, the IRAS 16293- 
2422 spectra do not show the wealth of spectral features 
due to other complex organics - although the confusion 
limit is clearly not yet reached. The limits on the abun- 
dances of molecules such as CH3OCH3 and HCOOCH3, 
derived if it is assumed that these molecules are located 
only in the warm inner part of the envelope, are higher 
than those found for high-mass protostars (Table |?]). Deep 
integrations down to ~5-10 mK are needed to verify the 
presence of these complex organic species. 

Typically, species such as SO2, OCS, HC3N and 
CH3CN are assumed to be 'second generation' products, in 
which case their abundance ratios can be used as 'chemical 
clocks' to determine the time since the evaporation of the 
ice mantles. For the density and temperature prevailing in 
the inner parts of the envelope around IRAS 16293-2422 
an age of ~10 4 yr is inferred from such chemical models, a 
value that is rather uncertain and sensitive to the cosmic 
ray ionization rate (e.g., Charnley 1997; FTatchell et al. 
1998). However, if the velocity field derived from the dust 
modelling for the Shu infall model is correct, the transit 
time for grains and molecules through the warm, dense 
region surrounding IRAS 16293-2422 is only several hun- 
dred years. This is not sufficient time for extensive second 
generation processing to result, and would naturally ex- 
plain the lack of large, complex organic species toward this 
source if their abundances are indeed found to be very low 
by subsequent deep searches. In such a scenario, the ob- 
served abundances would provide a unique opportunity to 
trace the chemical richness derived from previous stages of 
grain mantle and gas phase chemistry at sensitivities much 
higher than those achieved via infrared spectroscopy of icy 
grains. 



6.2. Alternative scenarios 

The above discussion demonstrates that the IRAS 16293- 
2422 data can be consistently interpreted with a 'hot core' 
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Fig. 9. Observed line widths (AT/) of different molecules 
as functions of their estimated rotational excitation tem- 
perature (T rot ). 



in which the ices are liberated from the grains by thermal 
evaporation above 90 K. However, one problem is that the 
size of the 'hot core' is only ~150AU in radius, compa- 
rable to the size(s) of the circumstellar disk(s), and much 
smaller than the binary separation. On this scale, depar- 
tures from spherical symmetry can be expected to play a 
role and the disk can shadow part of the inner envelope, 
keeping it much colder. Moreover, the IRAS 16293-2422 
envelope is known to be traversed by outflows. To what 
extent could these processes affect the interpretation? The 
data clearly indicate jumps in the abundances of 'first gen- 
eration' molecules known to be present in interstellar ices, 
but could these icy mantles also be removed by alternative 
mechanisms? 



Blake et al. (1994) and van Dishoeck et al. (1995 ) pro- 
posed that grain-grain collisions in the turbulent shear 
zones where the outflow interacts with the envelope can 
also be effective. This mechanism is observed for the class 
protostar L1157, where the outflow interaction can be 
spatially separated from the immediate protostellar en- 
vironment with single-dish telescopes. The inferred abun- 
dances by |Bachillcr fc Perez Gutierrez (1997 ) are included 
in Table [7] and are seen to also be close to the values found 
in the inner region of IRAS 16293-2422, especially for the 
sulfur-bearing species. 

An interpretation in which the ices are liberated by 
mild shocks or turbulence rather than thermal evapora- 
tion has two observational consequences. First, the lines 
of the 'first generation' ice mantle species are expected 
to be wider and have a different velocity structure than 
those of molecules located predominantly in the quiescent 
outer envelope. Figure ^ shows the observed line widths 
of different molecules as functions of their excitation tem- 



perature, where the values are taken from Blake et al 
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(1994) and [van Dishoeck et al. (1995| ). Molecules with 
clear 'jumps' in their abundances (SO2, CH3OH, CH3CN) 
have larger line widths and higher excitation tempera- 
tures than molecules which trace the outer envelope (CN, 
C 2 H, DNC, DCO+). The line widths observed towards 
IRAS 16293-2422 are significantly wider than what is typ- 
ically observed for class sources (<lkms _1 for C 18 
and C 17 0; J0rgensen et al. 2002) The observed widths 
up to Skms^ 1 could, however, possibly be explained by 
the infall model with thermal evaporation (Fig. |J) and 
are not clear-cut evidence for association with the molec- 
ular outflow, illustrating the difficulty in disentangling the 
contributions from various velocity components. Second, 
the spatial distribution of the 'first generation' molecules 
will be different. Molecules produced by thermal evapo- 
ration should be located within a — 150AU (1") radius, 
whereas in the shock scenario they are expected to coat 
the walls of the outflow(s). Such an 'X-type' geometry 
can extend over a much larger region, even though the 
total mass of warm gas may be comparable to that in the 
first model. Chemically, there is expected to be little dif- 
ference between the two scenarios, except perhaps in the 
'second generation' products if the liberation by shocks oc- 
curs in lower density or temperature gas. Observations at 
sub-arcsec resolution with the Smithsonian SubMillimeter 
Array (SMA) and the Atacama Large Millimeter Array 
(ALMA) are needed to distinguish these scenarios. 

Recently, Viti et al. (2001) have suggested, based on 
chemical modelling, that some molecular abundance ra- 
tios are strongly affected by the presence of a shock. In 
these models, the shocks not only liberate the ice man- 
tles but also drive high-temperature (— 2000K) reactions. 
In particular the combination of HCO/H2CO and NS/CS 
ratios may be suited for tracing the dynamical history of 
a hot core. For the HCO/H2CO ratio we derive an up- 
per limit of ~1, about an order of magnitude larger than 
any model predictions, whereas no NS data are available. 
Thus, the abundances derived from the data set cannot 
be used to constrain their origin in the inner hot region of 
IRAS 16293-2422. 

7. Conclusions 

The continuum emission radiated by dust grains in the 
circumstellar environment of recently born stars can be 
used, when supplemented by a detailed radiative transfer 
analysis, to derive the physical properties of such dusty 
envelopes. Observational constraints in the form of both 
the SED and resolved radial brightness distributions are 
required for a successful model. For the protostellar ob- 
ject IRAS 16293-2422 it is possible to model the observa- 
tional data using a single power-law distribution for the 
density structure. However, the molecular line profiles sug- 
gest that parts of the envelope are in a state of collapse. A 
dynamical age of ~l-3xl0 4 yr is derived with an annual 
mass accretion rate of roughly 4xl0 -5 M yr _1 . 

Once the physical structure of IRAS 16293-2422 is 
known, it is possible to determine its chemical proper- 



ties through detailed radiative transfer modelling of the 
observed molecular line emission. The abundances and 
quality of the fits for molecules like CO and CS further 
strengthen the adopted physical model. While the emis- 
sion from some molecules is well reproduced assuming a 
constant fractional abundance throughout the envelope, 
other species require a steep abundance gradient to be in- 
troduced at typically 90 K. The presence of such a jump 
for molecules like H2CO and CH3OH is interpreted as ev- 
idence of thermal evaporation of ices in the inner dense 
and hot regions of the envelope. Other molecules like 
HC3N, CH3CN and several sulfur-bearing molecules also 
require such drastic jumps in their abundance distribu- 
tions. In high mass protostars, these molecules are typi- 
cally cited as evidence of a rapid gas-phase chemistry ini- 
tiated by the evaporation of the ices. The 'hot core' region 
in IRAS 16293-2422 is, however, very small, only -150 AU 
in radius and comparable to the size of the circumstellar 
disk(s). This large change in physical size and the orga- 
nized velocity field surrounding low mass protostars lead 
to chemical time scales that only a small fraction of the 
dynamical age. Thus, it may be difficult for a full 'hot core' 
chemistry leading to complex species such as CH3OCH3 
to develop in class objects, even though the appropriate 
physical conditions are present. Alternative scenarios in 
which ices are liberated by grain-grain collisions in turbu- 
lent shear zones associated with the outflows need to be 
tested by higher angular resolution observations. 

Whatever their precise origin, the molecules located 
in the inner envelope of IRAS 16293-2422 can be incor- 
porated into the growing circumstellar disk(s) and be- 
come part of the material from which planetary bodies are 
formed. The molecular abundances derived here should 
provide an accurate reference point for comparison with 
the growing amount of data on protoplanetary disks and 
icy solar system objects such as comets. 
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